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Abstract . Tho dynamics of orbiting shallow flexible spherical ahell struc-* 
tures under the influence of control actuators la studied# Control lavs are 
developed to provide both attitude and shape control of the structure. It 
la seen that the elastic modal frequencies for the fundamental and lower 
modes are closely grouped due to the effect of the shell curvature. The 
shell is also assumed to be gravity stabilised by a spring-loaded dumbbell 
type damper attached at its apex. Control laws are developed based on the 
pole clustering technique and it is assumed chat tho dumbbell state informa- 
tion may not be directly observable. Numerical results verify that a signi- 
ficant savings in fuel consumption can be realized by using the hybrid shell- 
dumbbell system together with point actuators. Other results indicate that 
for the less robust systems Instability may result by not including the 
orbital and first order gravity-gradient effects in the plant prior to con- 
trol law design. 

Keywords . Modelling of orbiting flexible structures; pole placment; 
^deili^g errors I hybrid coucrol aystemn. 


INTRODUCTION 

Future proposed space missions would ‘involve 
large inherently flexible systems for use in 
communicatians , radlometry, and in electronic 
orbital based ma:?l systems. The use of very 
large shallow dish type structures to be em- 
ployed as receivers /reflectors for these 
missions has been suggested. In order to 
satisfy mission requirements control of the 
shape as well os the over-all orientation 
will be often required. The proposed paper 
is devoted to a study of the shape and ori- 
entation control of such an orbiting shallow 
spherical shell structure and, to the author^ 
knowledge, represents Che first such treat- 
ment of this subject. 

A related recent paper (Kumar and Bainum, 
1981) treated the dynamics and stability of 
a flexible spherical shell in orbit in the 
absence of active shape and orientation con- 
trol. For small amplitude elastic displace- 
ments and rigid rotational modal amplitudes, 
it was seen that the roll-yaw (out-of-plane) 
motions completely separate from the pitch 
(In-plane) and elastic motions. Furthermore, 
the pitch and only the axi-symmetric elastic 
modes are coupled within the linear range* 


*Research supported by NASA Grant NSG-1414, 
Suppl. 3. 


With the symmetry axis nondnally following 
the local vertical, the structure is gravi- 
tationally unstable due to an unfavorable 
moment of inertia distribution. A rigid 
dumbbell connected to the shell at its apex 
by a spring loaded double gimhallad joint 
with damping was proposed to gravitationally 
stabilize the structure (Fig. 1). It was 
noted that Che dumbbell motion could excite 
only those elastic modes having a single 
nodal diameter (Heissner, 1955) and that to 
completely damp the system transient motion 
in all of the Important lower frequency 
modes, the ixaa of an active control system 
would be required. 

The present paper represents an extension 
of the paper by Kumar and Bainxun (1981) to 
include in the mathematical model of Che 
dynamics Che effects of point actuators lo- 
cated at pre-selected pou,itions on the shell 
surface (Fig. 1). 


DEVELOPMENT OF MATHEMATICAL MODEL 
OF THE PLANT 

The mathematical model of an isotropic shal- 
low flexible spherical shell in orbit was 
developed by Kumar and Bainum (1981) under 
the assumption that the shell's elastic dis- 
placements were principally in the transverse 
direction (parallel to the symmetry axis) 
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Center of curvature 

Fig» 1. Shallow spherical isliell with duabbell and accuators 
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and wara amalX as aomparad with tha ochat 
charactarisclc dimanslons of the shalX. Tha 
assumption of shaXXownass fusrchar Insuras 
that tha ratio of tha dlspXacamant of tha 
ahaXX's apa^ point abova its base plana (U) 
is small as cooparad with tha radius In tha 
basa plana, a, (Fig* 1)* 


Tha resulting linaarlzad aquations of motion 
for tha rigid rotational and g^narlc elastic 
modas vara devalopad aa: 

♦"■wn,++a-fi,)i('' - c,/j u2 

Z Z Z S w 

0"-3ny9-2pj[ij^>Z/jy - Cy/jyU| 

c''+({S5“3)e„+26'lf“VM„Jt - 3 i5“Vm I 

n n n JL n l n 

+ (a-1,2,...- ) <1> 

vhara the variables and constant coafflclents 
are defined in the Appendix • (An 
order of magnitude analysis also indicated 
chat tha coupling batvaen the rigid (orbital) 
translationai nudes and thosa dascribad in 
Eq»(l) was extremely small for structures 
with characteristic lengths of 100m* so that 
these modes are essentially governed by the 
orbital machanico of the aystam mass center.) 

It was further ass timed that a dumbbeU could 
be attached by a spring loaded gimbal damper 
to the shell at its apex and could provide 
both gravitational stability of the uncon- 
trolled system as well as passive restoring 
and dissipative forces* 

The linearized equations of motion for the 
shallow spherical shelly-dumbbell system 
were developed ast 

r-nj-a+D^)^' - 
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Y'''<<y(l-tei)Y'+{3'H?ya-h;j,))Y+3(l+ny)8 

(I+C2 ) « (l+oj) }8+4a-nj) 4 

It is seen from £q* (1) that for the uncon- 
trolled system without the dumbbell that the 
out-of-plane roll-yaw motions are completely 
decoupled from the in-plane pitch (9) and 
elastic motions (Ca)* Within the linear 
range, only the sjdsymaetric elastic modes 
are coupled to the pitch motion. 
Futhermore from the analysis of Eq. (2) by 
Kumar and Bainum (1981) it was concluded 
that the dumbbell motion could excite only 
those elastic modes having a single nodal 
diameter t:nd that to completely damp the 
system transient motion in all of the impor- 
tant lower frequency modes, the use of an 
active control system would also be required 
However, it was hoped that a properly de- 
signed hybrid control system consisting of 
the passive dumbbell and active control 
actuators could provide satisfactory per- 
formance with a savings in fuel consumption 
as compared with the active thrusters oper- 
ating alone* 


The formulation of the uncontrolled dynamics 
assumes an a priori knowledge of the fre- 
quencies of ail the elastic imsdes to be 
incorporated within the system model. The 
frequencies (p) of the spherical shell are 
evaluated using the fallowing identities, 
as presented by Johnson and Reissner (1958) ; 


u “ [ ^ (p^p^)]^ (3) 


where the u's are calculated from 



for n “ 0,1 
(4) 


where n represents the number of nodal dia- 
meters (meridians), 2>»EhVl2(l-v2) , and p^ 
E/pR^. For n>l, Eq. (4) taust be replaced 
by a more complex form aa follows: 

F " FTIo' 

a 

Sjj(u)- 4n2(a2-l)(l-vHu[J„(ii)i;(w) 
]+(n+l) (1-v) [I^(u)- a 

-n^j^Cu) ]+W^Jjj(u) }{ (l-v)n 2 [vl^(w) 
-i^(y)]-w^l^(u)M(l-v)n2[uJ^(u) 
]+v^j;(u) Ka-v) ]ux^(ij) 
-n\(u)]-u\(u)} 
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and are Beaael functlorm of the flrac 

kind and^modlfiftd Beasel functions of the 
first kind, respectively, E<i. (4) or Eq. 

(5) is amcisfied by an infinite number of 
the paramacett for every value of n 
(j«l»2f •••«*)• For the sample calculations 
in this paper; we will consider only three 
such values of u (J”* It2f3) for the cases 
where n >* 0,1* 

The values of the natural frequencies and 
mode shape functions of the aadsymmetric 
toodes will be slightly modified by the 
presence of the dumbbells However^ f^r 
this applicationt an order of magnitude 
analysis for the system parameters involvedy 
indicates that the coupling between the 
axisymmetric modes and the rigid pitch mode 
is extremely weak and that> to a good first 
approximation f the small number of axi-* 
symmetric modes included can be considered 
independentJLy of all the rigid rotational 
modes, In view of this the axiaymmetric 
frequencies and mode shapes as given Johnson 
and Relssner (1958) are used here as a first 
approximation to the actual values in the 
presence of the dumbbell, the natural fra-* 
quencies and mode shapes of the other 
elastic modes characterized by nodal diama- 
ters (iMridians) remain unaffected by the 
presence of the dumbbell (Kumar and Bainumi 
1981). 

The point actuators are modelled as follows. 
An actuator located at (x^y^z) with com- 
ponents provides the following 

torques, ^ 

and the corresponding generic force in the 
nth mode, 

/ O-f* dm (7) 

n n 

For the shallow spherical shell it is - 
assumed that the major elastic displacement 
occurs In a direction norml to the base 
(y,z) plane. - l.e. ♦ Thus for a 

point actuator locatea at L3C,y,z) 

- 4“^<x,y.z)f^M^ W 

where if^*'\x,y,z) Is the a th modal shape 
function evaluated at (x,y,z), f^^ repre- 
sents the component of 7 in a direction 
normal to the base plane and is the mo- 
dal mass, here considered to be the total 
mass of the shell, since the shell mode 
shape functions have each been normalized 
with respect to the mass of the shell. 


NUMERICAL EXAMPLE 

As on example a large flexible shallow shell 
Is selected with the following dimensions 
(Fig, 1) and material properties: 


H n Im; a . lOOn; p " 27.68kg/m^; h - lea; 

V « 1/3; E - 0.744xl0^®kg/a2; R - 3000.0m 

Six elastic modes are included in the trun- 
cated model and are selected to have either 
no or a single nodal diameter# The number 
of nodal circles is varied from 1 to 3 such 
that the elastic modal frequencies are cal- 
culated using Eqs, (3) and (4), The results 
are summarized in Table 1# 

TABLE 1 - Elastic Modal Frequencies o f 
the SheU (Six Modes) 


of ttod^ * if of nodal P^(H^) 
diameters (n) circles (j) 


0 

1 

0 

1 

0 

1 


1 0,107593022 

1 0.107678417 

2 0.107946351 

2 0,108475664 

3 0,109516711 

3 0,111142560 


Xe is observed that all the frequencies are 
grouped with only a difference of 31 between 
maxiiitum and minimum values of the flexible 
frequencies considered in this model, this 
clwac grouping is due to the do m i n a n te of 
the curvature of the shell on the assumed 
model, and further emphasizes the importance 
of a careful consideration of potential mo- 
delling errors on the design of the control 
system. . 


The resulting state vectots for the system 
now take Che form: 

18x1 (9) . 

without the dumbbell; and 

22x1 (10) 

with the dumbbell. 

The equations of motion nsn be written in 
the state vector format; 

X « AX+BU (11) 

where the B matrix is evaluated using the 
location of the actuators, the values of the 
elastic modal functions at these locations 
(Johnson and Relssner (1958) and Itao and 
Crandall (1979)) and the direction cosines 
of the actuator thrust vectors. 


It is assumed that six actuators are posi- 
tioned along the surface of the shell (in 
the x ,2 plane as shown in Fig. 1), Force 
directions are selected so that when the 
actuators are operating torques will be pro- 
vided directly about each of the shell’s 
principal axes. The assumed actuator loca- 
tions and force directions are summarized 
in Table 2. 
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TABtE 2 - 

Actuator locations and Force 
Directiotui 

A cotq}arlson of the maximum force amplitude 
in each of the actuators and also of the 

Actuator Koi 


8 

X 

y 

t 

total force impulse required ia given in 
Table 3 for the two cases • It can be seen 

X) ^ 

90° 

0.28 


0 

0 

that although there la little difference in 
the maximum force amplitudas, there ia ap- 

2) 

90° 

0*57 

*2 

0 

0 

proximately a 25 percent savings la fuel 
consumption by using the shell-dumbbell sys' 

3) 

90° 

0*84 

h 

0 


tern* 

4) 

270° 

0*28 

h 

0 


TABIH 3 - Maximum Force Amplitudes of 
the Shell with 6 Actuators 

5) 

270° 

0.57 

*5 


0 

6) 

270° 

0*84 

^6 


Without With 

dumbbell dumbbell 


Tha locatloa of tho non«z«ro of 

tho ayaceoi niacrix; la llluatratad in Figa* 
2 and 3 for tha ayatam without tha duahbalX 
and containing tha dinobball* raapectivaly, A 
comparlaon of chasa figuras iULuatrataa the 
graatar coupling with tha rigid and flaadbla 
modaa that ia provided by tha dumbbell ♦ Nu-* 
oarical valuao for Che complaca A and B ma- 
tricae are listed in Bainua and co*-workers 
(1981) * It is seen that thw B matrix of tha 
shall with tha dumbbell haa the form; 

r 0 1 


Lb- J 

where 3* la the (9xd) lower portion of tha 
B matrix without the dumbbell* 



103.4UNt. 

103.410Nt. 

fj 

105.398 

105.398 


317*83 

317.83 

^4 

107.55 

107.55 


77*09 

79*09 


127.32 

127*318 

Total 

Force impulse; 

432*3 Mt-aec 

320*49 Ut-sec 
(25? savings 


with dumb-* 
bellO 


Control laws are selected using the pole 
\ clustering algorithm developed by Armstrong 
(1978) for the spherical shell without tha 
dumbbell* This aaioe control law is then 
applied to the case of the shell with tho 
dumbbell, i*e*- it is assumed that the dumb- 
bell position and rate information may not 
be directly observable and la not included 
in the control law* (Aa long as the dumb- 
bell-shell system is controllable and sta- 
bilirable, the dumbbell will return to its 
desired local vertical equilibrium orienta- 
tion after the transients have been removed* 
It is possible to design a control law for 
the du^bell-ehell system that results in a 
controllable-stabilisable system without 
including the dumbbell information within 
the control law*) 

A typical time history of the required 
control forces for the shell without the 
dumbbell is shown in Fig, 4 for non-dimen- 
sional Initial position displacements in all 
state components of 0«01» with the control 
law based on placement of all the poles so 
as to have a dimensionless negative real part 
of (-1*72)# This will provide a system cima 
constant of 461 seconds in all modes* A 
similar response was generated vsing the same 
control law for Che case of the shell with 
the dumbbell where the dumbbell inertia 
ratioa (cj^^c,) were aasi^dj:o be^ 0*9, Che 
dumbbell spring constants, 
and dacping coefficients selected^ so as to 
provide 0*1 critical damping. 


Uony current investigationa of the shape 
and orientation control of flexible orbiting 
structures do not incorporate the effects of 
the gravity-gradient and orbital dynamic 
coupling into the linear model of the plant* 
-i.e. the poles of the rigid rotational 
modes of the open loop system are at the 
origin* A study was made to determine the 
effects of omitting these terms in the devel- 
opment of control laws* Control laws based 
on pole placement ware first developed based 
on the shell model which did not Include the 
gravity-gradient and orbital coupling terms 
in Che plant modal* The control laws thus 
developed were then inserted into the pre- 
viously developed models which contain both 
first order gravity-gradient and orbital 
dynamic coupling terms* Typical results are 
illustrated In Figs* 5 and 6* In all cases 
studied, there is a general tendency for some 
of the poles of the rigid modes to ahi^c 
towards the imaginary axis when the gravity- 
gradient and orbital effects are superimposed 
Into the plant* For the less robust systems 
(Fig* 5) instability may result* In general, 
there is no noticeable shift in the poles 
corresponding to the flexible modes* As 
expected, for ttkft more robust systems, the 
relative effect of this shift is less appar- 
ent (Fig* 6), but at the expense of greater 
control force effort (Table 4)* 
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Fig, 2* Location of the noa-tero aleoaata of the 
matrix of shell without the dumbbell. 
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Fig. 3. Location of non-zero elements of the 
matrix of the shell with the dumbbell. 
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O (a) dttsiigatd cloiAd loop poles 

model does not include orbital 
find grA\rit 7 ** 5 radisa& effects 




Fig* 6, Shift in closed loop poles due to orbital and gra'/ity-gradient effects-designed 
response time - 400 secs. 
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TABLE 4 - roro Attplitu4<i« 

Control law dftvoloptd without orbital and 
9r4vity**grAdi<mt nffactA in tha plane and 
than oppUad to a taodal including th«i«» 



Without 
orbital 1 
gravity- 
gradient 

With orbital 6 

gravlty- 

gr.di.nt 

*1 

516.40Nt. 

565.63Ht. 

fj 

73.30 

164.76 

*3 

239.50 

321.11 


117.54 

U7.07 

^5 

132.45 

163.84 


146.82 

431.61 

Total 

fore. 

lopuls.: 

7020.8Nc-sec 

19373.0Kt-g.c 


CONaUSIONS 

OriAneacion and ahape control of an orbiting 
ahalXov spharical ahaXX syatom may bo acconr* 
pliahod by uaing appropriatoly positionod 
actuators on tho surfaco of tho shall* A 
gimballad spriag-*loadad dumbboll dampor 
connocted at tho shall Vs apox con provide 
gravitational stabilization together with 
a source of passive damping. A significant 
savings in fuel consumption can be realized 
by using the combined active end passive 
control systems. For less robust systems 
instability may result by not including 
the orbital and gravity-gradient effects 
in the plant prior to the design of the 
active control law. 
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APPENDIX - NOMENCLATURE 


a 

C 


C,Cy,C^ 


: modal lusplitudes 
I base radius of shell 
: external torques with components 

3 coefficients of viscous damping 


-(n) f,(a) . 


m. 


(n) 


ay ^ ( 0 . 0 ). 

r.sp.ctlv«ly 


jT 


( 0 . 0 ). 


-(mn) „ (mn). . y , c(m)-(o) 




D 

E 


a 

(n) 



c(n»)c(n) 
z z 


‘ 

respectively 

: flexural rigidity 

: modulus of elasticity 

: modal component of external 

forces 


vol 


Xc*x“^dv 


moment of inertia of the dumb- 
bcU 


J ,J 
y* z 


principal momenta of Inertia 
of the undeformad shell 


k.ky.St, 


: torsional spring constants 

* respectively 

; characteristic length (-a, the 
base radius) 

: modal mass of n^^ mode 


radius of curvature of the 
shell 
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i 

coordliuet of difftraticial 
artA oa tb« «urf^i abovt 
tbM bane plana 

t 

t 

tijxm 

0 

1 

polar angltt of particular 
location on ahall 


• 

» 

anglft EMtautad in baaa plana 
of oball 

Y.0 

: 

duobball daflacclon anglta 

‘a 

: 

Kj% 


t 

pitch 1 yav and roll anglaa» 
raapactivaly 

V 

t 

Poisson *9 ratio 

p 

9 

mass density 

T 

t 

(ii^t (dimansionltss tioa) 

u 

: 

body angular va/lf.Alty vactori 
or 

“c 

t 

orbit angular valoclty 

< 

9 

th 

natural fre<iuaacy of n node 


t 

r*ap«ctlv«l 7 

(•)•( )' 

S 

3t ' ^ ' rtap«ct±v«ly 

.M 

• 

Cr«n*v*rfla conponanc at th« 
n ^ modal sbapa function 
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